Omega from the skewness of the cosmic velocity divergence

We propose a method for measuring the cosmological density parameter $\Omega$ from the statistics of the divergence field, $\theta \equiv H^{-1} \div v$, the divergence of peculiar velocity, expressed in units of the Hubble constant, $H \equiv 100 h km/s/Mpc$. The velocity field is spatially smoothed over $\sim 10 h^{-1} Mpc$ to remove strongly nonlinear effects. Assuming weakly-nonlinear gravitational evolution from Gaussian initial fluctuations, and using second-order perturbative analysis, we show that \propto -\Omega^{-0.6} ^2. The constant of proportionality depends on the smoothing window. For a top-hat of radius R and volume-weighted smoothing, this constant is $26/7-\gamma$, where $\gamma=-d\log / d\log R$. If the power spectrum is a power law, $P(k)\propto k^n$, then $\gamma=3+n$. A Gaussian window yields similar results. The resulting method for measuring $\Omega$ is independent of any assumed biasing relation between galaxies and mass. The method has been successfully tested with numerical simulations. A preliminary application to real data, provided by the POTENT recovery procedure from observed velocities favors $\Omega \sim 1$. However, because of an uncertain sampling error, this result should be treated as an assessment of the feasibility of our method rather than a definitive measurement of $\Omega$.Comment: 16 pages + 2 figures, uuencoded postscript file, also available by anonymous ftp from ftp.cita.utoronto.ca in directory /cita/francis/div_skewness, CITA 94-1

Three-Point Statistics from a New Perspective

Multipole expansion of spatial three-point statistics is introduced as a tool for investigating and displaying configuration dependence. The novel parametrization renders the relation between bi-spectrum and three-point correlation function especially transparent as a set of two-dimensional Hankel transforms. It is expected on theoretical grounds, that three-point statistics can be described accurately with only a few multipoles. In particular, we show that in the weakly non-linear regime, the multipoles of the reduced bispectrum, $Q_l$, are significant only up to quadrupole. Moreover, the non-linear bias in the weakly non-linear regime only affects the monopole order of these statistics. As a consequence, a simple, novel set of estimators can be constructed to constrain galaxy bias. In addition, the quadrupole to dipole ratio is independent of the bias, thus it becomes a novel diagnostic of the underlying theoretical assumptions: weakly non-linear gravity and perturbative local bias. To illustrate the use of our approach, we present predictions based on both power law, and CDM models. We show that the presently favoured SDSS-WMAP concordance model displays strong ``baryon bumps'' in the $Q_l$'s. Finally, we sketch out three practical techniques estimate these novel quantities: they amount to new, and for the first time edge corrected, estimators for the bispectrum.Comment: 5 pages 6 figures, ApL accepte

Extended Perturbation Theory for the Local Density Distribution Function

Perturbation theory makes it possible to calculate the probability distribution function (PDF) of the large scale density field in the small variance limit. For top hat smoothing and scale-free Gaussian initial fluctuations, the result depends only on the linear variance, sigma_linear, and its logarithmic derivative with respect to the filtering scale -(n_linear+3)=dlog sigma_linear^2/dlog L (Bernardeau 1994). In this paper, we measure the PDF and its low-order moments in scale-free simulations evolved well into the nonlinear regime and compare the results with the above predictions, assuming that the spectral index and the variance are adjustable parameters, n_eff and sigma_eff=sigma, where sigma is the true, nonlinear variance. With these additional degrees of freedom, results from perturbation theory provide a good fit of the PDFs, even in the highly nonlinear regime. The value of n_eff is of course equal to n_linear when sigma << 1, and it decreases with increasing sigma. A nearly flat plateau is reached when sigma >> 1. In this regime, the difference between n_eff and n_linear increases when n_linear decreases. For initial power-spectra with n_linear=-2,-1,0,+1, we find n_eff ~ -9,-3,-1,-0.5 when sigma^2 ~ 100.Comment: 13 pages, 6 figures, Latex (MN format), submitted to MNRA

Self-similarity and scaling behavior of scale-free gravitational clustering

We measure the scaling properties of the probability distribution of the smoothed density field in $N$-body simulations of expanding universes with scale-free initial power-spectra, with particular attention to the predictions of the stable clustering hypothesis. We concentrate our analysis on the ratios $S_Q(\ell)\equiv {\bar \xi}_Q/{\bar \xi}_2^{Q-1}$, $Q \leq 5$, where ${\bar \xi}_Q$ is the averaged $Q$-body correlation function over a cell of radius $\ell$. The behavior of the higher order correlations is studied through that of the void probability distribution function. As functions of ${\bar \xi}_2$, the quantities $S_Q$, $3 \leq Q \leq 5$, exhibit two plateaus separated by a smooth transition around ${\bar \xi}_2 \sim 1$. In the weakly nonlinear regime, {\bar \xi}_2 \la 1, the results are in reasonable agreement with the predictions of perturbation theory. In the nonlinear regime, ${\bar \xi}_2 > 1$, the function $S_Q({\bar \xi}_2)$ is larger than in the weakly nonlinear regime, and increasingly so with $-n$. It is well-fitted by the expression $S_Q= ({\bar \xi}_2/100)^{0.045(Q-2)}\ {\widetilde S}_Q$for all$n. This weak dependence on scale proves {\em a small, but significant departure from the stable clustering predictions} at least for n=0$and$n=+1$. The analysis of$P_0$confirms that the expected scale-invariance of the functions$S_Q$is not exactly attained in the part of the nonlinear regime we probe, except possibly for$n=-2$and marginally for$n=-1$. In these two cases, our measurements are not accurate enough to be discriminant.Comment: 31 pages, postscript file, figure 1 missing. Postscript file including figure 1 available at ftp://ftp-astro-theory.fnal.gov:/pub/Publications/Pub-95-256-

Galactic dust polarized emission at high latitudes and CMB polarization

We estimate the dust polarized emission in our galaxy at high galactic latitudes, which is the dominant foreground for measuring CMB polarization using the high frequency instrument (HFI) aboard Planck surveyor. We compare it with the level of CMB polarization and conclude that, for angular scales $\le 1^{\circ}$, the scalar-induced CMB polarization and temperature-polarization cross-correlation are much larger than the foreground level at $\nu \simeq 100 GHz$. The tensor-induced signals seem to be at best comparable to the foreground level.}Comment: Latex document, 4 pages, 2 figures, to appear in "Fundamental parameters in Cosmology", Rencontres de Moriond, 199

Error estimation for the MAP experiment

We report here the first full sky component separation and CMB power spectrum estimation using a Wiener filtering technique on simulated data from the upcoming MAP experiment, set to launch in early 2001. The simulations included contributions from the three dominant astrophysical components expected in the five MAP spectral bands, namely CMB radiation, Galactic dust, and synchrotron emission. We assumed a simple homogeneous and isotropic white noise model and performed our analysis up to a spherical harmonic multipole lmax=512 on the fraction of the sky defined by b>20 degrees. We find that the reconstruction errors are reasonably well fitted by a Gaussian with an rms of 24 $\mu$K, but with significant deviations in the tails. Our results further support the predictions on the resulting CMB power spectrum of a previous estimate by Bouchet and Gispert (1999), which entailed a number of assumptions this work removes.Comment: 5 pages, 3 color figures, version accepted in A&A Letter

Simulations of the Microwave Sky and of its ``Observations''

Here follows a preliminary report on the construction of fake millimeter and sub-millimeter skies, as observed by virtual instruments, e.g. the COBRA/SAMBA mission, using theoretical modeling and data extrapolations. Our goal is to create maps as realistic as possible of the relevant physical contributions which may contribute to the detected signals. This astrophysical modeling is followed by simulations of the measurement process itself by a given instrumental configuration. This will enable a precise determination of what can and cannot be achieved with a particular experimental configuration, and provide a feedback on how to improve the overall design. It is a key step on the way to define procedures for the separation of the different physical processes in the future observed maps. Note that this tool will also prove useful in preparing and analyzing current (\eg\ balloon borne) Microwave Background experiments. Keywords: Cosmology -- Microwave Background Anisotropies.Comment: 6 pages of uuencoded compressed postscript (1.2 Mb uncompressed), to appear in the proceedings of the meeting "Far Infrared and Sub-millimeter Space Missions in the Next Decade'', Paris, France, Eds. M. Sauvage, Space Science Revie